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Determination of Leak and Respiratory Airflow 

Field of the Invention 

5 The invention relates to methods and apparatus for the determLnauon of Jeakage 

airflow and true respiratory airflow, particularly during mechanical ventilation. 

The airflow determination can be for a subject who is either spontaneously ar_ 
riqn-spontaneously breathing, or moves between these breathing states. The invention is 
10 especially suitable for, but not limited to, normally conscious and spontaneously 
breathing human subjects requiring long term ventilatory assistance, particuiarly during 
sleep. 

]5ackground of the Invention 

IS 

In this specification any reference to a "mask* is to be understood as including 
all forms of devices for passing breathable gas to a person's airway, including nose 
masks, noise and mouth masks, nasal prongs/pillows and cndocracheai or tracheostomy 
cubes. 
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During mechanical ventilation, breathable gas is supplied for example via- a 
mask, at a pressure which is higher during inspiration and lower during expiration, h 
is useful to measure the subject's respiratory airflow during mechanical ventilation to 
assess adequacy of trcatmeni, or to control the operation of the ventilator. . 
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Respiratory airflow is commonly measured with a pne\imotachograph placed in 
the gas delivery path between the mask and the ventilator. Leaks. between the mask and 
the subject arc unavoidable. The pneumotachograph measures d\e sum of the 
respiratory airflow plus the tlow through die leak. If Uie' instantaneous flow through die 
30 leak IS known, die respiratory airflow can be calculated by subtracting the flow through 
the leak from die flow ai the pneumotach. 

Known methods to correct for the flow through the leak assume (i) that tile 
leak is substantially constant, and (ii) that over a sufficiendy long time, inspiratory and 
35 expiratory respiratory airflow will cancel. If these assumptions arc met, the average 
How dirough the pneumotach over a sufficiently long period will equal the magniwdc 
of the leak, and the true respiratory airflow may dien be calculated as described. 



The knowa method is only correct if the pressure at the mask is conscani. If 
the., mask pressure varies with time (for example, in the case of a vemilaror), 
assumption (i) above will be invalid, and the calculated respiratory airflow will 
therefore be incorrect. Tliis is shown markedly in Figs, la- If. 

Fig. la shows a trace of measured mask pressure in bi-level CPAP treatment 
between about 4 cm H2O on expiration and 12 cm H2O on inspiration. Fig. lb shows 
a trace of true respiratory airflow in synchronism with the mask pressures. At time 
21 seconds a mask leak occurs, resulting in a leakage flow from the leak tliat is a 
function of die treatment pressure, as shown in Fig. Ic. Tne measured mask flow 
shown in Fig. Id now includes an offset due to the leak flow. Tlie prior art method 
then determines the calculated leak flow over a number of breailis, as shown in Fig. le. 
The. resulting calculated respiratory flow, as the measured flow minus die calculating 
leak flow is shown in Fig. If, having returned to the correct mean value, however is 
incorrecdy scaled in magnitude, giving a false indication of peak positive and negative 
airflow. 

Another prior art arrangement Is disclosed in European Publication 
No. 0 714 670 A2, including a calculation of a pressure-dependent leak component. 
The methodology relies on knowing precisely the occurrence of the start of an 
inspiratory event and the start of the next inspiratory event. In Q±cr words, the leak 
calculation is formed as an average over a known breadi and applied 10 a subsequent 
breath. 

This method cannot be used if the moment of start and end of the previous 
breaus are unknown. In general, it can be difficult to accurately calculate the time of 
start of a breath. This is panicularly the case immediately following a sudden change 
in the leak. 

Furthermore, the method will not work in tlic case of a subject who b; making 
no respiratory efforts, and is momentarily not being ventilated ai all, for example 
during an apnea, because for the duration of the apnea tliere is no start or end of brcatli 
over which to make a calculation. 

The present invention seeks to provide a determinaiion of leak flow and irue 
respiratory airflow, accounting for the variations in flow through a leak as a funciion of 
pressure. 



Summary of the Invesdon 



The invention discloses a mediod f©r determining instantaneous Icalc flow at a 
mask having a leak path during mechanical ventilation, the mcrhod comprising the siep.s 
of: . 

(a) determining inscauianeous airflow at the mask; 

(b) determining instantaneous pressure at the mask; 

(c) estimating non-linear conductance of said leak path as the low-pass 
filtered instantaneous airflow divided by the low-pass filtered square root of the 
instantaneous pressure; and 

(d) determining said instantaneous leak flow lo be said conductance 
multiplied by the square root of the said instantaneous pressure. 

The invention, further discloses a method for determining instanraneous 
respiratory airflow for a subject receiving breathable gas by a mask and in the presence 
of any masJc leak, the method comprising the steps of: 

(a) determining instantaneous airflow at the mask: 

(b) determining instantaneous pressure at the mask: 

(c) estimating non-linear conductance of said leak path as the low pass 
filtered instantaneous airflow divided by tlie low pass filtered square root of tlie 
instantaneous pressure; 

(d) determining instantaneous leak flow to be said conductance multiplied 
by die square root of die said instantaneous pressure; and 

(e) calculating the respiratory airflow as the instantaneous airflow minas 
the instantaneous leak flow. 

The invention yet further di'scloses apparatus for determining respiratory 
airflow for a subject rcceivmg breathable gas by a mask and in die presence of any 
mask leak, the apparatus comprising: 

transducer means located at or proximate die mask to determine instantaneous 
mask airflow and pressure; and 

processing means for estimating non-linear conductance of said leak patii as the 
low pass filtered instantaneous airflow divided by the low pass filtered square root of 
die instantaneous pressure, deteirmining instantaneous leak flow to be said conductance 
multiplied by the square root -of the said insianumeous pressure, and calculating ihc 
respiratory airflow as the instantaneous airflow minus the instantaneous leak flow. 

The invention yet further discloses apparatus for providing continuous positive 
airway pressure treatment or mechanical ventilation, the apparatus cofnprising: 



a lurbine for the generation of a supply of breathable gas; 
a gas delivery tube having connection with the turbine; 
a mask having, connection to the delivery vubc to supply said breatliable gas to 
a subject's airway; 

transducer means located at or proximate tlic mask to determine instantaneous 
mask airflow and pressure; 

processor means for estimating non-linear conductance of said leak path as the 
low pass filtered instantaneous airflow divided by the low pass filtered square root of 
the instantaneous pressure, determining instantaneous leak flow to be said conductance 
m.ultiplied by the square root of the said instantaneous pressure, and calculating tlie 
respiratory airflow as the instantaneous airflow minus the instantaneous leak flow; and 

control means to control ihe flow generator to, in turn, control the mask 
pressure and/or mask airflow on the basis of the calculated respiratory airflow. 

The invention yet furtlier discloses a computer program for. executing the steps 
referred to above. 

In one preferred form, time constants of tlie low pass filtering are dynamically 
adjusted dependent upon sudden changes in the instantaneous leak tlow. 

Embodiments of the invention provide advantages over the prior art. There is 
no need to know when transitions between respiratory phases occurs. The 
independence from knowledge of the subject's respiratory state has the important result 
diat the leak flow calculation is accurate in apneic (i.e. no flow) instances on die part of 
the subject or the mechanical ventilator. 

Brief Description of the DrairtDgs 

Embodiments of the invention will now be described with reference to the 
accompanying drawings in which: 

Figs, la-lf shows traces of pressure and airflow from which respiratory 
airflow is calculated in accordance with a prior art method: 

Figs. 2a and b show schematic diagrams of iwo embodiments of ventilatory 

assistance apparatus; 

Fig. 3 is a block flow diagram of a metliod Ibr determining instantaneous 

respiratory airflow; and. 

Fig. 4 shows traces of pressure, airflow and odier variables From which 

respiratory airflow is calculated; 



Fig. 5 shows a schematic diagram of ventilatory assistance apparatus of another 
embodiment; 

Fig. 6 sIjows a fuzzy membership function for the calculatioii of the' extent Aj 
to which the time t^x since die most recent positive going zero crossing of the calculated 
respiratory airflow is longer than die expected time Tj; 

Fig. 7 shows a fazzy membership function for the calculation of die cxtenr B| 
to ^'hich the calculated inspiratory respiratory airflow fRg^p is large positive: 

Fig. 8 shows a fuzzy membership function for the calculation of the extent Ag 
to which the time t^E since the most recent negative going zero crossing in tfie 
calculated respiratory airflow is longer tlian the expected time Tg; 

Fig. 9 shows a fuzzy membership funaion for the calculation of the extent B£ 
to which die respiratory airflow f^gsp ^ '^g^ negative; and 

Fig. 10 shows die relation between an index J and time constant t. 

Detailed Description of Preferred Embodiments 

Fig. 2a shows mechanical ventilation apparatus 10 embodying tlie invention. 

The subject/patient wears a nose mask 12 of any known type. The subject 
equally could wear a ftice mask or nasal prongs/pillows, or alternatively have an 
endotracheal tube or tracheostomy tube in place, A turbine/blower 14, operated by a 
mechanically coupled electrical motor 16, receives air or breadiable gas at an inlet 18 
thereof, and supplies the breadiable gas at a delivery pressure to a delivery tube/hose 20 
having connection at the odier end thereof widi die nose mask 12. Breaduble gas thus 
is provided to the subject's, airway for the purpose of providing assisted respiration, 
witli the subject's expired breath passing to atmosphere by an exhaust 22 in die delivery 
cube 20, typically located proximate to the mask 12. 

A pneumotachograph 24 is placed in tlic deliver^' tube 20 between the mask 12 
and the exhaust 22 to provide two pressure signals, P2 and P(, across the 
pneumotachograph, each passed by hoses 28,30 to a differential pressure sensor 32. A 
determination of the flow of gas in the mask 12 is made the differencial pressure, P2- 
?!, resulting in a flow signal fj. The mask pressure, Pi. also is passed to a pressure 
sensor 34 by a tapped line 36 taken from rhe respective hose 28, to generate a delivery 
pressure signal, p^^, output from die pressure sensor 34. 

Bodi die flow signal, f^j, and die pressure signal, p^, are passed to a 
microcontroller 38 where they are sampled for subsequent signal processing, typically 
''at a race of 50 Hk. 



The microcontroller 38 is prograxnmed to process the flow aiid pressure signals 
<fd» ^m) '° produce an output conffol signal, yo. provided to an electronic motor servo- 
controUer 42 that, in turn, produces a motor speed control output signal, Vq. This 
signal is provided to the motor 16 to control the rotational speed of the turbine 14 and 
provide the desired treatment pressure, ?2, at tlie nose mask 12. 

The motor servo-controUer 42 employs a negative feedback control technique . 
that compares the actual del Wery pressure, in ±e form of die signal p^. with Ihe 
control signal y©. For convenience, this control stratagem may be independent of 
operation of the microcontroller 38. 

Operaiion of the controUiiig of the microcontroller 38. so far as a calculation 
of respiratory airflow is concerned, broadly is as follows. In a sampled manner, the 
conductance of any "mask leak is calculated, then the instantaneous flow through the 
leak is calculated. The flow through die leak is subtracted from the total mask flow to 
calculate the true instantaneous respiratory airflow. 

Fig. 2b shows an alternative embodiment of a system for determining true 
respiratory ^airflow during mechanical vejitiiation. The mechanical vcntUation system 
10* of Fig. lb differs from chat of Fig. la firstly in that the microcontroller 38 plays no 
part in control of die ventilator 50. rather only receives and data processes the 
electrically transduced mask pressure and flow signals Pm- fd » dcxBrmme and generate 
the instantaneous respiratory flow fRgsP- The ventUator 50 has an internal drive signal 
provided by an oscillator 44. The motor servo controller also may or may not receive 
the mask pressure signal pn, as a form of feedback control. Indeed, the ventilator 50 
can be realised by any convenient form of known generic ventilation device. 

The controlling software resident within the microcontroller 38 performs ihe 
following steps in determining the respiratory airflow as broadly described above, as 
also shown in the flow diagram of Fig. 3. 

The word "average" is used herein in die most general sense of the result of a 
low pass filtering step, and is not confined to an aridvnetic mean. 

1. Repeatedly sample the mask airflow to give a sampled signal fMASK. and the 
mask pressure to give a sampled signal PmaSK-. for sample at intervals of T=2D 
milliseconds. (Steps 50.52). 
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2. Calculate the average leak, LP(L), as being the resulc of low pass filtering the " 
airflow, fMASK' "^'^^ ^ ^nie constant of 10 seconds. (Step 54), 

3. Calculate tlie average of the square root of the mask pres.'^ure, LPC^/Pj^^ask). as being 
5 the result of low pass filtering the square root of the mask pressure, PmaSK» a 

time constant of 10 seconds. (Step 56). 

4. Calculate the conductance, G, of any leak (Step 58). from the equation: 

10 G =LP(L)/LP(V7mask) 

5. Calculate the. instantaneous leak flow. fLEAK, through the leak (Step 60), from the 
equation: 

If there is ao leak flow, tlie value of LP(L) will be equal io zero, as will G and hence 
fLEAK. Thus the methodology is valid also where lealc is equal tn zero - no leak. 

10 At this juncture the leak flow has been determined, such as would be desired for a leak 
flow detector. If desired, the instantaneous respiratory airflow can be subsequently 
determined by the following step. 

6. Calculate the instantaneous resphatory airflow, fRESP.. by subtracting the 
2s instantaneous leak firora die mask flow (Step 62): 

^ESP = fMASK- fLEAK 

Figs. 4a-4h illustrate the methodology of the embodiment described above with 
reference to Fig. 2b. At time, t = 21 sec, a continuing leak of approximately I I/sec is 

30 inrroduced. Fig. 4e shows the mean mask flow. Fig. 4f represents the calculated 
conductance G, from which the mask leak flow can be estimated as shown in Fig, 4g, 
Finally, Fig. 4h shows how the calculated respiratory airflow recovers within 
approximately 30 seconds, and, importandy, gives the correctly scaled (true) niagnitude. 
of airflow. 

35 . • 

Widi regard to setting the instantaneous output signal yQ, ilie microcontrolter 
broadly executes the following steps: 



r 
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7. If ±e calculated true respiratory airflow fRJESP is above a threshold, for example 
0.05 L/sec, yo is set to a value corresponding to an inspiratory pressure, PtnSP- 
Otherwise yo is set co a value corresponding :o an expiratory pressure, PexP- 
general, PiNSP '^'gl^c^ -han PexP* "se of continuous positive airways 

.pressure, PeXP ^INSP- (Step 66). 

It is to be understood that many odicr methods of determining yo from fMASK 
may be used in step 7, for example as described in the text Principles and Practice of 
Mechanical Yemilalion. edited by Martin J. Tofain (McGraw Hill Inc, 1994). 

In order to control ventilation, it is necessary to measure the subject's 
ventilation. In the presence of a leak, the ventilation delivered by the assisted 
ventilation apparatus is greater than the ventilation delivered to the subject. Known 
devices which servo-control ventilation cope with this by collecting the exhaled air 
stream with a complex system of valves, and then measuring the exhaled ventilation. 
This is inappropriate for devices for use in a domestic setting during sleep, because of 
the attendant weight, complexity, and expense. The embodiment described compensates 
for the leak by continuously measuring the nonlinear conductance of the leak, and 
allowing for the instantaneous flow dirough the leak as a function of pressure. 

Fig. 5 shows an alternate arrangement for ventilatory assistance apparatus 10' 
embodying the invention. In this arrangement, the pneumotachograph 24' is interposed 
between the lurbine 14 and tlie delivery hose 20. 

This arrangement removes the pressure sensing hoses and pneumotachograph 
from tlie region of die mask 12. The pressure at the mask, PmaSK- calculated from 
die delivery pressure at die turbine 14, and from the pressure drop down the air 
delivery hose 20, which for any panicular delivery hose is a known function of the flow 
at the pneumotachograph 24. Further, die microconcroUer 38 must also calculate die 
flow dirough die mask from die flow at the mrbine 14 less die flow tiirough die exhaust 
22, which for any particular exhaust is a known function of die pressure at the mask 12. 

In more detail, diis involves the steps of, FirsUy measuring the pressure P3 at 
the turbine 14 with die pressure sensor 34 to produce an electrical signal pj. Next die 
differential pressure p4-p3 is measured across the pneumotachograph 24* by die 
differential pressure sensor 32 to produce an electrical signal ff a sumpled manner, 
pt and. ft are digitized to yield die sampled njrbine pressure and flow signals 
PtURBINE and FturbINE- 



Tlie pressure at the mask PmaSK *® sampled airflow at the mask fjvfASK 
12 are calculated from the turbine pressure PtURBINE ouUet of the 

turbine F7URBINE ^ foHows: 

1. Calculate the pressure drop APjube down the air delivery tube 20, from the 
flow at the outlet of the turbine FtURBINE- 

^PtUBE = sign (FturBINE) ^ K] (FtuRBINE)^ + ^2 PTURBfNE . 

where and K2 are empirically determined constants, and sign (x) is I for 
x>Q and -I otherwise. 

2. Calculate the pressure at the mask, PmASK' ^ pressure at the turbine 
p-pURBINE the pressure drop APjuBE ^^^^ delivery tube 20: 

Pmask = Pturbine - ^PtUBE' 

3. Calculate the flow". fsxHAUST' thrpugii exhaust 22, from the pressure at 
the mask PmaSK- 

^EXHAUST = sign(PMASK) ^ ^3 >/ahs(PMASK) 

where K3 is determined eii^)irically, 

4. Calculate the flow, fMASK> ^"^ ^-^^^ turbine 14 less 
the flow through the exhaust 22: 

fMASK = ^TURBINE " ^EXHAUST- 

The foregoing embodiments describe low-pass filtering of both die 
instantaneous airflow and the square root of the instantaneous pressure with a time 
constant t of 10 seconds. This time constant, t, can be advantageously dynamically 
adjustable. 

If the conductance of the lealc suddenly changes, then Uie calculated 
conductance will initially be incorrect, and will gradually approach ihc correct value at 
a rate which will be slow if the time constant of the low pass filters is long, and fast if 
the time coa;tant is short. Conversely, if .die impedance of the leak is steady, the 
longer the time constant the more accurate die calculation of. the instantaneous leak. 



- 10- 



Therefore, it is desirable to lengthen the time constant if it is certain that the leak is 
steady, reduce the time constant if it is certain that the leak has suddenly changed, and 
to use imermediately longer or shoner time consanis if it is intermediately certain tiiat 
die leak is steady. 

If there is a large and sudden increase in the conductance of the leak, tlien the 
calculated respiratory airflow will be incorrect. In particular during apparent 
inspiration, die calculated respiratory airflow wUl be large positive for a time that is 
large compared with the expected duration of a normal inspiration. Conversely, if thefe 
is a sudden decrease in conductance of the leak, then during apparent expiration the 
calculated respiratory airflow will be large negative for a time that is large compared 
with the duration of normal expiration. 

Therefore, an index of the degree of certainty that the leak has suddenly 
changed is derived, such that the longer the airflow has been away from zero, and by a 
larger amount, the larger the index; and the time constant for die low pass, filters is 
adjusted to vary inversely with the index. In operation, if tliere is a sudden and large 
change in the leak, the index wiU be large, and the time constant for the calculation of 
the conductance of the leak will be small, . allowing rapid convergence on the new value 
of the leakage conductance. Conversely, if tlie leak is steady for a long time, the index 
will be small, and the time constant for calculation of the leakage conductance will be 
large! enabling accurate calculation of the instantaneous respiratory airflow. In the 
spectnim of intermediate situations, where tlie calculated instantaneous respiratory 
airflow is larger and for longer periods, the index will be progressively larger, and the 
time constant for the calculation of the leak will progressively reduce. For example, at 
a moment in time where it is uncertain whether the leak is in fact constant, and the 
subject merely commenced a large sigh, or whether in. fact diere has been a sudden 
increase in the leak, the index will be of an intermediate value, and die time constani 
for calculation of the impedance of the leak will also be of an intermediate value. One 
advantage is that some corrective action will occur very early. 

Anodier advantage is diat there is never a moment where the leak correction 
algorithm is "out of control" and needs lo be restarted, a.s described for prior art 
European Patent Publication No. 0 714 670 A2. 

In a preferred embodiment, the above index is derived using fuzzy logic. -Hie 
fuzzy extent A, to which the airflow has been posuivc for longer than expected is 
calculated from the time i^X since the last positive-going zero crossing of the calculated 
respiratory airflow signal, and die expected duration T^of a normal inspiration for the 



particular subject, using the fuzzy membenhip function shown in Fig. 6. The fuzzy 
extent Bj to which the airflow is large and positive is calculated From the instantajieous 
respiratory airflow using the fuzzy membership function shown in Fig. 7. The 
instantaneous index of the degree of certainty tiiat the leak has suddenly increased is 
calculated as the fuzzy- intersection 0«ser) of Aj and Bi. 

Comparable calculations are performed for expiration as follows. The fuzzy 
extent A£ to which the airflow has been negative for longer than expected is calculated 
from the time t^s since the last negative-going zero crossing of the calculated 
respiratory airflow signal, and Tg, the expected duration of a typical expiration for the 
particular subject, using the membership function shown in Fig. 8. The fuzzy e^xsm 
Bg to which the airflow is large negative is calculated from the inscantaneous 
respiratory airflow using the fuzzy membership function shown in Fig. -9. The 
instantaneous index 1^ of the degree of certainty that tiae leak has suddenly decrea.sed is 
calculated as the fuzzy intersection of Ag and 

The instantaneous index I of the extent to which dierc has been a sudden 
change in die leak (either an increase or a decrease) is calculated as the fuzzy union 
(larger) of indices Ij and Ig. The instantaneous index I is Chen passed through a peak 
detector followed by a low pass filter with a time consiant of. for example 2 seconds, to 
yield, the desired index J. Thus if index I becomes momentarily large, index J will be 
initially large and remain so for a few seconds. The time constant t for ihe low pass 
filters used in the calculation of the conductance of the leak is then adjusted to vary 
inversely witli the index J, as shown in Fig. 10. For example, if the expected duration 
of a normal respiratory cycle were 4 seconds the time consiant is set co 10 seconds if 
the index J is zero, (corresponding to complete certainty that the leak Is steady), and to 
1 second if the index J is unity (corresponding to complete certainty that the leak is 
suddenly changing), and to intermediate values for intermediate cases. 

The embodiments described refer to apparatus for the provision of veniiiatory 
assistance, however, it is to be understood diat tlie invention is applicable co all forms 
of meclianical ventilation and apparata"; for the provision of continuous positive airway 
pressure treatment. Tlie apparatxis can be for die provision of a constant treatment 
pressure, multi-level (IPAP and EPAP) treatment or autosetting (adjusting) rreatmem or 
other forms of mechanical ventilation, including Proportional Assist Ventilation (PAV) 
as caught by M Younes in die above-noted text. 
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The methodology described can be implemented in the form of a computer 
program that is executed by the microcontroller described, or by discrete combinational 
logic elements, or by analog hardware. 



